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ABSTRACT 

The calculation of the functions is presented in a way that leads to the 
calculation of the stiffness matrix of an element with linear reinforced concrete 
brackets. The numerical development is oriented to three cases, namely: i) Beam 
with a carton considered as a single element; ii) Portfolio on the left with constant 
section; iii) Constant section and cartouche on the right. In this way the Structural 
Planner has two options to perform static or seismic analysis; the first without 
considering knot in the middle of the beam; and the second considering knot in the 
middle of the beam. 

Calculation tables are then presented to determine the elements of the 
flexibility and stiffness matrix of a beam with two linear brackets, for the indicated 
cases. Finally, the new functions of the CEINCI-LAB Computer system for the static 
and dynamic analysis of reinforced concrete structures with cartéis are indicated. 
Key words: Repair and reinforcement of structures. 

Keywords: Reinforced Concrete Linear Posters. CEINCI-LAB Computer system 


ESTRUCTURAS CON CARTELAS LINEALES DE 
HORMIGÓN ARMADO Y NUEVAS FUNCIONES DE 

CEINCI-LAB 

RESUMEN 

Se presenta el cálculo de las funciones de forma que conducen al cálculo 
de la matriz de rigidez de un elemento con carteias lineales de hormigón armado. El 
desarrollo numérico está orientado a tres casos, a saber: i) Viga con carteia 
considerada como un solo elemento; ii) Carteia a la izquierda con sección constante; 
iii) Sección constante y carteia a la derecha. De esta manera el Proyectista 
Estructural tiene dos opciones para realizar el análisis estático o sísmico; la primera 
sin considerar nudo en la mitad de la viga; y la segunda considerando nudo en la 
mitad de la viga. 

Luego se presentan tablas de cálculo para determinar los elementos de la 
matriz de flexibilidad y de rigidez de una viga con dos carteias lineales, para los 
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casos indicados. Finalmente se indican las nuevas funciones dei sistema de 
computación CEINCI-LAB para el análisis estático y dinâmico de estructuras de 
hormigón armado con carteias. 

Palabras Claves: Carteias Lineales de Hormigón Armado. Sistema de computación 

CEINCI-LAB 


1. INTRODUCCIÓN 

La tendencia es construir edificios con luces grandes, con parqueaderos 
subterrâneos en los que se pueda circular sin temor atoparse con alguna columna 
y tener todas las instalaciones: eléctricas, electrónicas, mecânicas, agua contra 
incêndios, etc., en forma visible para poder dar mantenimiento en el futuro. 

Para lograr satisfacer estos requerimientos se deben tener vigas con 
carteias como la indicada en la figura 1. En el análisis estructural se requiere mayor 
sección en los extremos con respecto al centro de luz, por ese motivo las carteias 
es una solución óptima y se deja la parte central para el paso de todas las 
instalaciones. 

b b 



Lv 


Lv 


1^_L_j 

Figura 1 Viga de hormigón armado, con carteias lineales simétricas y nomenclatura 
utilizada para definir la geometria de la viga. 

Sea b el ancho constante de la viga; h la altura en el centro de luz; H la 
altura en el extremo de la viga; L v la longitud de la carteia; L la longitud total de la 
viga, ver figura 1. 

En este artículo se deducen las funciones de forma y se determinan los 
términos de la matriz de flexibilidad y de rigidez, para el modelo i ndicado en la figura 
1, destacando que con este modelo no se puede considerar nudo en la mitad de la 
viga, sino que se lo debe analizar como un solo elemento a la viga. 

Ahora, si se desea considerar nudo en la mitad de la viga, se debe 
determinar las funciones de forma para los dos casos que se indican en la figura 2, 
a la izquierda se presenta el modelo de una viga con carteia al lado izquierdo y la 
gráfica adjunta es para el caso de carteia a la derecha. En la figura 2, se indica 
también la nomenclatura utilizada, solo difiere en la longitud L. 
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Tipo = 0 Tipo = 1 

Figura 2 Modelo que permite considerar nudo en la mitad de la viga; carteia a la 
izquierda y carteia a la derecha. 


Para los tres casos: i) Viga con dos carteias lineales simétricas (Figura 1); 
ii) Viga con carteia a la izquierda; iii) Ciga con carteia a la derecha (Figura 2) se 
determina en primer lugar las funciones de forma para el sistema 1, que no 
considera desplazamientos como cuerpo rígido, para el caso de carteia izquierda y 
lado recto este sistema se indica en la figura 3. El nudo inicial tiene un apoyo fijo y 
el nudo final un apoyo deslizante (rodillo). 

h 

í 

b 




Figura 3 Viga con carteia a la izquierda en el sistema de coordenadas 1 , apoyo fijo 
a la izquierda y rodillo a la derecha. 

En Aguiar (2020), Lamar y Fortoul (2007) se presenta el desarrollo numérico 
para obtener la matriz de rigidez de un elemento de sección constante, en 
coordenadas locales (sistema 2) y en coordenadas globales (sistema 3) a partir de 
la matriz de rigidez en el sistema 1. Es importante tener presente esto debido a que 
únicamente se muestra en este artículo el cálculo en primer lugar de las funciones 
de forma para una viga en sistema 1, luego se obtiene la matriz de flexibilidad y la 
inversa de esta matriz es la de rigidez. 

Una vez que se tiene la matriz de rigidez en el sistema 1 (apoyado-rodillo) 
el problema se considera resuelto, ya que la matriz de paso dei sistema 1 ai 2 
denominada 7\_ 2 y la matriz de rotación T 2 _ 3 es la misma que se tiene para un 
elemento de sección constante. Aguiar (2020). 
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2. VIGA SIN NUDO INTERMÉDIO 

Para encontrar la función de forma 0 x (x) en el sistema 1, se aplica un 
momento unitário en el nudo inicial, ver figura 4, consecuencia de esto se tienen las 
fuerzas de corte en los apoyos 1/L. 



'M 


Figura 4 Momento aplicado en el nudo inicial de vigacon carteia lineal simétrica. 
Sea: 


a =H 




h = a H 

• Para: 0 < x < Lv 


K= PL 


-Lv 



x = 0; y — H 
x = Lv, y = h 


m = 


h-H 


y — H = m(x — 0) 
H-h 

y = H - 7 -x 


/w = 4(„ 
,w = n( 


Lv 


H-h 


Lv 


H-aH 
H -——x 


pL 


bH 3 i 1 — a ' 3 

'“'tfI 1 — w x 

bH 3 

Ia = 


12 


Luego, la ecuación que define la variación de la inércia es la siguiente: 
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/ 1 — a \ 

IM = la{ 1 - —x) 


fSL 


• Para: Lv < x < L — Lv 


Io = 


b * h á 


Io = 


12 

b * (a//) 3 


12 

Io = a 3 Ia 


• Para: L-Lv <x <L 



-I v— 


x = L — Lv, y = h 
x = L; y = H 


m = 


H-h 

Lv 


y — H = m(x — L) 

y = H - — (L - x) 
Lv v J 


I(x) = — 
12 L 


H 


H -aH 


,3 


Lv 


<L-x) 


I(x) = 


bH 


3 


1 — a 


,3 


I(x) = Ia 


12 /?L 

1 — a 


■ (L - x) 


PL 


(L - x) 


La ecuación de momentos para cuando se tiene un momento unitário en ei 
nudo inicial es la siguiente. 


M(x) = X -~ 1 

El objetivo es calcular en primer lugar la ordenada de la elástica v(x), que 
es la función de forma 0 x (x) y después el giro 0(x) que al evaluarlo en x = 0, reporta 
el valor de Ai y al calculado en x = L, se tiene el valor de /l2> que son elementos 
de la matriz de flexibilidad para el sistema 1. (Aguiar, 2020). Para el efecto se tiene 
la ecuación de resistência de materiales que se indica a continuación. 
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v(x ) 


d 2 v M (x) 
dx 2 EI(x ) 


Para: 0 < x < (3L 


d z v 

dx 2 


r-i 


d 2 v 


E * la ( 1 ~ L pr x ) 

1 x-L 


dx 2 E * Ia * L 


d 2 v jS 3 * L 2 


x — L 


dx 2 E * Ia (J3L — (1 — a)x) 3 
Resolviendo por integraciones sucesivas se tiene que: 


dv (L 3 * L 2 


dx E * Ia 


2 (a — l)x + (/? — a + 1 )L 
2(a — 1 ) 2 ((a — l)x + /?L) 2 


+ A 


( 2 ) 


/? 3 * L 2 
E * Ia 


((2a — 2)x + 2/?i) ln((a — l)x + /?i) + (/? + a — 1 )L 


2 (a — l) 3 ((a — l)x + /?L) 


+ Ax + B 


Para: f!L < x < L - (IL 


d 2 v 


x z-' 


dx 2 E * a 3 * Ia 
d 2 v 1 


dx 2 E * a 3 * Ia * L 


(x-L) 


(3) 


Aplicando integraciones sucesivas se tiene que: 


dv 1 /x 2 \ 

— =---- Lx + C 

dx E * a 6 * Ia * L \ 2 J 


(4) 


1 (x 3 x 2 

v(x) = ---- L -b Cx + D 

E * a 3 * Ia * L \ 6 2 


(5) 
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• Para: L - pL < x < L 


d 2 v 
dx 2 


- 1 


E * Ia 


1 — a 

~w~ 


(L - x) 


d 2 v p 3 * L 2 x — L 

dx 2 E * Ia (J 3 L — (1 — cr)(L — x)) 3 


Integrando por integraciones sucesivas setieneque: 


dv p 3 * L 2 


2 (a - l)x- (P + 2 a- 2 )L 
dx E * Ia \2 (a — l) 2 ((/? + a — 1 )L — (a — l)x) 2 

p 3 * L 2 /ln((l — a)x + (P + a — 1 )L) 


+ E ' 


( 6 ) 


y(x) = 


E * Ia 


(a — l ) 3 
/?L 


2 (a — l) 3 ((cr — l)x — (/? + cr — 1)L) 


+ E* x + F 


( 7 ) 


Para encontrar las constantes de integración se tienen las condiciones de 
frontera que se incidan a continuación. 


x — 0; v = 0 


B = 

x = L; v = 0 

[ln((l — cr)L + (/? + cr — 1)L) 


(2 pi) ln (PD + (p + a - 1)L _ 
° “ 2(a — l) 3 (/?L) 

(2 pi) ln(/?L) + (p + a - 1)L 


2 (cr — l) 3 (/?L) 


( 8 ) 


0 = 


/?L 


(cr— l) 3 2(cr — l) 3 ((cr—1)L — (/? + cr — 1)L) 

ln(/?L) /?L 


+ E*L + F 


0 = 


(cr-1) 3 2 (a — 1) 3 (—/?L) 


+ E L+F 


0 = 


ln(/?L) 


+ ■ 


F = 


(cr — l) 3 2(cr —l) 3 
ln(/?L) 1 


(cr — l) 3 2(cr —l) 3 


+ F*L + F 


- EL 


( 9 ) 
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x = (3L;vi = vd 


p 3 *L 2 ( ((2a — 2)/ZL + 2/ZL) ln((a — 1)/ZL + /ZL) + (J3 + a — 1)L 
E * Ia * \ 2(a - l) 3 ((a - 1 )/ZL + JT) 

+ apl + b\ 


p 3 *L 2 

ui =-* 

E * Ia 


p 3 * Zr 

ui =-* 

E * Ia 


(2a/ZL) ln(a/ZL) + (/Z + a — 1)L 
2(a - l) 3 ((a - lj/ZL + JE) 

(2a/Z) ln(a/ZL) + (/? + a — 1) 


+ A/ZL + B 


2 (a — l) 3 ((a — 1)/? + /Z) 


+ A/ZL + B 


vd = 


E * a 3 * Ia * L 


(/ZL) 3 (/ZL) 2 

—- L — -+ C/ZL + D 


E * a 3 * Ia * L 


(/ZL) 3 (/ZL) 2 

^t^-l^-^ + c/zl + zz 


* L 2 
E * Ia 


(2a/Z) ln(a/ZL) + (/Z + a — 1) 
2 (a — l) 3 ((a— 1)/Z + /Z) 


+ A/ZL + B 


(/ZL) 3 (/ZL) 2 

^7^-l^-- + c/zl + zz 


a" 


P 3 * L 3 


(2a/Z) ln(a/ZL) + (/Z + a - 1) \ 

2 (a — l) 3 ((a — 1)/Z + /Z) + P ) 


D = (a/ZL) 3 


(2a/Z) ln(a/ZL) + (/Z + a — 1) \ 

2(a _ 1)3(a/?) + A/ZL + ZZ J 


(/ZL) 3 _ r (/ZL) 2 
6 + 2 


- C£L 


( 10 ) 


x = (iL)Qi = Gd 


/Z 3 * L 2 

2(a - 1)/ZL + (/Z - a + 1)L 

E * Ia* 

2(a - l) 2 ((a - 1)/ZL+/ZL) 2 

/Z 3 *L 2 

■ 2a/ZL -2 PL + PL - aL + L 

' E * Ia* 

2(a - 1) 2 (a/ZL) 2 


+ A 


+ A 


8i = 


/Z 3 *L 2 
E * Ia 


— 2 a/ZL + /3L + aL — L 
2 (a — 1 ) 2 (a/ZL) 2 
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p 3 * L 2 

9i = --* 

E * Ia 

Gd = 


-2 a/3L + /3L + aL — L 


2 (a - 1 ) 2 (a/?L) 2 

r (PL) 2 


+ 4 


E * a 3 * Ia * L \ 2 


PL 2 + C 


p 3 * L 2 
E * Ia 


-2 a/3L + (IL + aL — L 


2 (a - l) 2 (a/?L) 2 
—2 a/3L + /3L + aL — L 


+ 4 




2 (a — l) 2 (a/?L) 2 


+ 4 


E * a 2 * Ia* L \ 2 

1 /(/?D 2 


- PL 2 + C 


PL 2 + C 


1 A P L ) 7 

A = - -— —-/?L 2 + C + 

(a/?L) 3 \ 2 P ' 


(a/?L) 3 V 2 

2a/?L - pL-aL + L 
2 (a — l) 2 (a/?L) 2 


x = L — pL, vi = vd 

v(x) 

I 

ui = 


3 2 

X X 


E * a 3 * Ia* L 


, ,-L-h Cx + D 

E*a i *la*L\6 2 


(L - /?L ) 3 (L - PL ) 2 
--——-L--í-^-+ C(L -pL) + D 


vi = 


L 5 -2pL 5 +2p 2 L 3 -p 3 L 3 L 3 -2pL 6 + p 2 L 


r3 I n2 j3 


E * a 3 * Ia * L 

— PL) + D 
1 

ui = 


+ C(L 


E * a 3 * Ia * L 


2PL 3 - p 3 L 3 -2L 3 


+ C(L - pL) + D 


P 3 * L 2 

vd =-* 

E * Ia 


ln((l — a)(L — pL) + (J3 + a — 1 )L) 


(a - l) 3 

_ PL _ 

2(a - l) 3 ((a - 1 )(L - pL) - (/? + a - 1)L) 

+ F 


+ £(L - PL) 


vd 


P 3 *L 2 

E * Ia 


ln (a/?L) 

(a — l) 3 + 


1 

2(a — l) 3 (a) 


+ £(L — PL) + F 
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2PL 3 - p 3 L 3 -2L 3 


+ C(L - pL) + D 


(apLy 


ln (a/?i) 


(a — l) 3 2 (a — l) 3 (a) 


+ e{l - pi) + F 


( 12 ) 


x = L — pL\ 9i = 9d 
8i = 

6i = 


E * a 3 * Ia * L 


(L - pLY 


- L(L - pL) + C 


E * a 3 * Ia * L 


L 2 - 2/3L 2 + p 2 L 2 


L(L - /?L) + C 


9i = 


p 2 ! 2 - L 2 


+ C 


E * a 3 * Ia * L\ 2 

2 (a - 1 )(L - pi) - (J3 + 2a - 2 )L 
E * Ia \2(a — 1) 2 ((/? + a — 1)L — (a — 1)(L — /?L)) 


p 3 *L 2 

8d = --* 


-+ E* 


P 3 *L 2 


8d = 


*( 


-2a PL + PL 


E * Ia \2(a — l) 2 (a/?L)‘ 


+ E* 


a 3 * L 


1 (tl^ +c ) = tiEj 

V 2 / E * Ia V 


-2apL + /?L 
E * Ia V2(a — 1 ) 2 (a/?L)‘ 


+ £"* 


e . = A-Í ^ 2 ~ l2 + c )+ 


(a/?L) 3 


2 (a — 1 ) 2 (a/?L) 2 


(13) 


Para hallar el valor de las constantes, se despeja todas en función de C. 
• Para encontrar el valor de la variable D, se trabaja con la ecuación (10) 


D = (jotpE) 3 


(2aP)\n(apL) + (p + a - 1) \ 

2( a _ 1 )s ( a p^ +ApL + Bj 


m 3 , m 2 

6 +L 2 


-CPL 


Se reemplaza el valor de B (ecuación 8) 


D = (cr/?L) 3 


( (2a/?)ln(a/?L) + (/? + a — 1) aoj 

\ 2(a-map) +AljL 

(2pL) ln(/?L) + (J3 + a — 1)L\ (/?L) 3 (pL) 2 

2(a — l) 3 (pL) ) 6 +L 2 


Se sustituye el valor de A (ecuación 11) 
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D = (a/?L) 3 


D = (afiL) 3 


ln(a/?L) (/? + a - 1) ln(/?L) (/? + a - 1) 

(a - l) 3 2(a - l) 3 (a£) + (a - 1)3 + 2(a - 1) 3 (/?) + ** 

QSL) 3 (jSL) 2 

ln(a) (/? + a — 1) 2a — 1 

~F \ n 4" 


(a — l) 3 2(a —l) 2 (a/?) 2(a—l) 2 a 2(a — l) 2 (a) 2 /? 

w 2 


D = (afiL) 3 


Q3L) 3 (/?L) 2 

- —fo— + L —~Y - CfiL + (3L 

(a/3 -Va 2 — 2) 


+ ■ 


2a 


2 

ln(a) 


-/?Z/ +C 

m 3 


2(a —l)(a) 2 /? (a — l) 2 a (a—l) 3 

, (/?£) 2 
+ PL ‘ H 


+ L 


m 2 


■- /?Z/ + C )-CpL 


D = (a/?L) 3 


(a/? + a 2 — 2) 2a — 1 ln(a) 
2(a — l)(a) 2 /? + (a — l) 2 a (a — l) 3 
ÍPL) 2 


m 3 ,, m 2 


+ (3L 


~PL‘ 


Nótese en la ecuación anterior que se elimina la variable C, con lo que se 
obtiene el valor de D, directamente. 


D = ( a/3L ) 3 


(a/? + a 2 — 2) 2a — 1 


+ - 


ln(a) 


m 3 ,, m 2 , (iSi ) 3 

- I Lj T r-, 


D = ( a/SL ) 


2(a—l)(a) 2 /? (a —l) 2 a (a —l) 3 

ln(a) \ (/?L) 3 /? Z L 


- /? 2 L 3 


(a/? + a 2 — 2) 2a — 1 
2(a — l)(a) 2 /? (a — l) 2 a (a — l) 3 


2,3 


+ - 


(14) 


• Para F se reemplaza la ecuación 13 en la ecuación 9: 


ln(/?L) 1 L (p 2 L 2 -L 2 ir \ L(2a/?L — /?L) 

(a — l) 3 2(a — l) 3 (a/SLyy 2 + C J 2(a - l) 2 (a/?L) 2 

InQSL) 1 L (/3 2 L 2 -L 2 ir \ 2a-1 

(a — l) 3 2(a — l) 3 (a/?L) 3 ^ 2 + C ) 2(a - l) 2 (a) 2 /? 


Utilizando la siguiente relación se obtiene la constante C: 
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3/?L 3 ~p 3 L 3 - 2L 3 
—--+ C(L — /?L) + D 



i 


ln(a/?Z,) 1 

(a — l) 3 2 (a — l) 3 (a) 


+ E*(l - fll) + F 


3/3L 3 - p 3 L 3 - 2L 3 

—---+ C(L — /?L) + D 

6 

(a/?L) 3 

~~ í 

ln(a/?L) 1 ln(a) 1 

(a — l) 3 2(a — l) 3 (a) (a — l) 3 2(a —l) 2 (a) 

/?L íp 2 L 2 -L 2 \ 2a-1 

“ (a/?L) 3 V 2 + “ 2(a — l) 2 (a) 2 


3/?L 3 - p 3 L 3 - 2L 


+ C(L — /?L) + D 


(apL)' 


ln(a) 


2a — 1 


PL(P 2 L 2 -L 


(a — l) 3 2(a — l) 2 (a) 2(a — l) 2 (a) 2 


+ C 


3/?L 3 - p 3 L 3 - 2L 


+ C(L — /?L) + D 


( a/3L )' 


ln(a) 


2a — 1 


P 3 L 3 - pL 3 


[(a — l) 3 2(a — l) 2 (a) 2(a — l) 2 (a) 2 J 


+ C/?L 


3 PL Ó - p 3 L 3 - 2L 


+ C(L — pV) + D + CPL 


(iapL )' 


ln(a) 


2a — 1 


(a — l) 3 2(a — l) 2 (a) 2(a — l) 2 (a) 2 


P 3 L 3 -PL 3 


[D + CL] = 


0 apL ) 3 

1 


ln(a) 1 2a — 1 

(a — l) 3 2(a — l) 2 (a) 2(a — l) 2 (a) 2 


P 3 L 3 - L 3 


3 


(a/?L) 3 

2 ln(a) 

6 ap +a 2 -2p-2 

i 3P 2 L 3 + 2L 2 - \p 3 L 3 
~h 

L 

(a - l) 3 

2 (a — l)(a) 2 /? 

6 
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• Teniendo C. Las constantes de integración se resumen en: 
1 


A = 


(J 3 L ) 2 , \ 2 aBL - BL-aL + L 

-/ 3 L 2 + c) + 


(a( 3 L ) 3 \ 2 J 2 (a — l) 2 (a/?L) 2 

(2 pi) ln(/?L) + 0? + a - 1)L 


B = 


2(a — l) 3 (J 3 L) 


C = 


(«Aí-) 3 

2ln(a) 

6a/? + a 2 — 2/? — 2 

* 

L 

í(«-l) 3 

2 (a — 1 )(a) 2 /? 


3/? 2 L 2 + 2L 2 - 4/? 3 L 


3 j 3 


+ ■ 


n , m ^3 , («/? + a: 2 — 2) 2a - 1 ln(a) \ (/?L) 3 /? 2 L 3 

D = (a/5L) 3 * ( —- — + 7 - 7T7 - 7 - 777 I + ■ 


2(a —l)(a) 2 /? (a —l) 2 a (a — l ) 3 


r = -L-(g!t!^ +c 1+ 


F = - 


(ctpL) 3 

ln(/?L) 1 


2 (a — 1 ) 2 (a/?L) 2 


L //? 2 L 2 - L 2 


(a — l) 3 2(a — l) 3 (a/?i) 3 


+ C - 


2a — 1 


2(a — 1 ) 2 (a) 2 /? 


Una vez que se han determinado las constantes de integración, se tiene 
que las funciones de forma son: 

• Para: O < x < Lv 

P 3 *l 2 


010 ) = 


E * Ia 


((2a — 2)x + 2/?i) ln((a — l)x + /?i) + (/? + a — 1)L 


2 (a — l) 3 ((a — l)x + /?L) 


+ Ax + B 

• Para: Lv < x < L — Lv 

010 ) = 


3 2 

X X 


, . — L -b Cx + D 

E * a 6 * la* L \ 6 2 


• Para: L — Lv < x < L 


I 3 3 *L 2 (\n{(l-a)x+ (/? +a- 1)L) 

01(x) = ——7— * 


E * Ia 


(a — l) 3 
/?L 


2 (a — l) 3 ((a — l)x — (/? + a — 1)L) 


+ E x + F 


A partir de los resultados encontrados, en Aguiar (2020) se deducen las 
funciones para ei caso de sección constante, como comprobación de lo realizado. 
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En base a los resultados obtenidos, se han elaborado, las tablas 1 y 2 que 
permiten el cálculo de los elementos de la matriz de flexibilidad y de rigidez, 
respectivamente en el sistema 1. 

Los elementos de la parte superior de la tabla 1, son los coeficientes que 
sirven para el cálculo dei valor de Ai ~ A 2 (por ser carteias simétricas) y los 
términos de la parte inferior (valores negativos) sirven para calcular fl2 — Al- 


Tabla 1 Coeficientes a* y /?*, para obtener la matriz de flexibilidad en vigas con 
carteias lineales simétricas. Fuente: Bonilla (2019) _ 


p\a 

0,95 

0,9 

0,85 

0,8 

0,75 

0,7 

0,65 

0,6 

0,55 

0,5 

0,45 

0,4 

0,5 

0,306 

0,279 

0,254 


0,209 

0,186 

0,165 

0,146 

0,128 

0,111 

0,095 

wm 

-0,157 

-0,148 

-0,139 


-0,121 

-0,112 

-0,103 

-0,094 

-0,085 

-0,077 

-0,068 


0,45 

0,308 

0,284 

0,260 

0,238 

0,217 

0,197 

0,178 

0,160 

0,144 

0,128 

0,113 

RI 

-0,159 

-0,151 

-0,143 

-0,136 

-0,128 

-0,121 

-0,113 

-0,106 

-0,098 

-0,091 

-0,084 

■S9 

0,4 

0,310 

0,288 

0,267 

0,247 

0,227 

0,209 

0,192 

0,175 

0,160 

0,145 

0,132 


-0,160 

-0,154 

-0,148 

-0,141 

-0,135 

-0,129 

-0,123 

-0,117 

-0,111 

-0,105 

-0,099 


0,35 

0,313 

0,293 

0,274 

0,255 

0,238 

0,221 

0,206 

0,191 

0,177 

0,164 

0,151 


-0,162 

-0,157 

-0,152 

-0,147 

-0,142 

-0,137 

-0,132 

-0,127 

-0,122 

-0,117 

-0,113 


0,3 

0,315 

0,297 

0,281 

0,265 

0,249 

0,235 

0,221 

0,207 

0,195 

0,183 

0,172 


-0,163 

-0,159 

-0,155 

-0,151 

-0,148 

-0,144 

-0,140 

-0,137 

-0,133 

-0,129 

-0,126 


0,25 

0,318 

0,303 

0,288 

0,274 

0,261 

0,248 

0,236 

0,225 

0,214 

0,204 

0,194 

0,185 

-0,164 

-0,161 

-0,159 

-0,156 

-0,153 

-0,150 

-0,148 

-0,145 

-0,142 

-0,140 

-0,137 

-0,135 

0,2 

0,320 

0,308 

0,296 

0,284 

0,274 

0,263 

0,253 

0,244 

0,235 

0,226 

0,218 

0,211 

-0,165 

-0,163 

-0,161 

-0,160 

-0,158 

-0,156 

-0,154 

-0,152 

-0,151 

-0,149 

-0,147 

-0,145 

0,15 

0,323 

0,314 

0,304 

0,295 

0,287 

0,279 

0,271 

0,264 

0,257 

0,250 

0,244 

0,238 

-0,166 

-0,165 

-0,164 

-0,163 

-0,161 

-0,160 

-0,159 

-0,158 

-0,157 

-0,156 

-0,155 

-0,154 


Ai - A 2 ~ a * e . 


fi2 - Ai - £ * E , 


(15) 


Al sacar la inversa de la matriz de flexibilidad, se tiene la matriz de rigidez, 
para el sistema de coordenadas 1 y se obtienen los coeficientes: k* (valores 
superiores de tabla 2) y a* (valores inferiores de tabla 2) 


Ai = k 22 = k* E -± 


(16) 


k 12 — Ai — a 


E In 
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Tabla 2 Coeficientes para encontrar la matriz de rigidez en vigas con carteias 
lineales. 


p\a 

0,95 

0,9 

0,85 

0,8 

0,75 

0,7 

0,65 

0,6 

0,55 

0,5 

0,45 

0,4 

0,5 

4,447 

4,975 

5,603 

6,359 

7,169 

8,423 

9,858 

11,695 

14,103 

17,341 

21,838 

28,335 

2,288 

2,635 

3,059 

3,582 

4,136 

5,062 

6,128 

7,529 

9,411 

12,007 

15,708 

21,192 

0,45 

4,424 

4,920 

5,507 

6,207 

7,053 

8,088 

9,374 

10,999 

13,096 

15,866 

19,633 

24,939 

2,282 

2,620 

3,031 

3,534 

4,157 

4,941 

5,941 

7,240 

8,960 

11,294 

14,552 

19,257 

0,4 

4,397 

4,859 

5,399 

6,036 

6,797 

7,715 

8,838 

10,232 

11,993 

14,261 

17,252 

21,304 

2,272 

2,596 

2,986 

3,459 

4,038 

4,757 

5,658 

6,807 

8,296 

10,261 

12,914 

16,587 

0,35 

4,367 

4,788 

5,275 

5,843 

6,509 

7,299 

8,246 

9,394 

10,804 

12,560 

14,780 

17,638 

2,258 

2,562 

2,923 

3,355 

3,875 

4,508 

5,285 

6,250 

7,462 

9,004 

10,995 

13,606 

0,3 

4,332 

4,708 

5,136 

5,625 

6,190 

6,844 

7,609 

8,509 

9,577 

10,855 

12,394 

14,264 

2,239 

2,517 

2,841 

3,221 

3,670 

4,202 

4,838 

5,602 

6,527 

7,655 

9,037 

10,743 

m 

4,292 

4,616 

4,979 

5,385 

5,843 

6,360 

PI 

7,614 

8,376 

9,248 

10,248 

11,394 

2,215 

2,460 

2,740 

3,060 

3,428 

3,853 


4,912 

5,572 

6,339 

7,232 

8,269 

0,2 

4,246 

4,513 

4,806 

5,126 

5,477 

5,862 

6,284 

6,747 

7,255 

7,812 

8,420 

9,081 

2,185 

2,390 

2,619 

2,874 

3,158 

3,475 

3,828 

4,222 

4,659 

5,145 

5,682 

6,272 

0,15 

4,194 

4,400 

4,619 

4,852 

5,100 

5,364 

5,643 

5,938 

6,250 

6,577 

6,919 

7,276 

2,148 

2,309 

2,482 

2,669 

2,870 

3,087 

3,320 

3,568 

3,834 

4,115 

4,413 

4,725 


3. CARTELAS PARA VIGAS CON NUDOS INTERMÉDIOS 


Con este modelo no se puede resolver estructuras en las cuales se 
considera un nudo en la mitad de la viga, por lo que en este apartado se presentan 
nuevos formulários para encontrar la matriz de rigidez en carteias de la forma 
indicada en la figura 2; ei modelo de la izquierda es para cuando la carteia se halla 
a la izquierda (Identificada como Tipo=0) y ei modelo de la derecha para cuando la 
carteia se halla a la derecha (Tipo=1). 

3.1 Cálculo de funciones de forma para carteia al lado izquierdo 

En la figura 5, se presenta la viga con carteia a la izquierda en el sistema 1 
(apoyo fijo-apoyo móvil) con su nomenclatura y a la derecha se muestra la viga en 
la que se aplica un momento unitário en el nudo inicial, consecuencia de este 
momento son las reacciones que valen 1/L; pero ahora L es la mitad de la longitud 
de la viga. De tal manera que se ha calcular la función de forma 0 í (x) en forma 
similar al desarrollado en el apartado anterior pero ahora es más fácil. 


Por didáctica se repiten las siguientes ecuaciones: 

h= a H L v = (3 L 
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Figura 5 Cálculo de las funciones de forma para carteia al lado izquierdo para 
momento aplicado en nudo inicial. 


• Para: 0 < x < Lv 



x = 0; y — H 
x = Lv, y = h 


m 


■ 

Lv 


y — H = m(x — 0) 


y = H 
b 


H-h 


Lv 


-Lv- 


' w = n( H ~ 

,M =à( 


H-h 


Lv 

b i H-aH ' 3 
t- 


(LL 


bH 3 ( 1 — a \ 3 

'^TTV—fTV 

bH 3 

Ia = - 

12 

l 1 — a X 
;w = la (l - — x) 


• Para: Lv < x < L 


Io = 


Io = 


(3L 

b * h 3 
12 

b * (a//) 3 


12 

Io = a 3 Ia 


La ecuación de momentos para cuando se tiene un momento unitário en ei 
nudo inicial es la siguiente. 

MÇx) =l~ 1 

d 2 v M (x ) 
dx 2 EI(x ) 
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• Para: 0 < x < pL 


d 2 v 
dx 2 


- 1 


d 2 v 


E * Ia 
1 


(it) 


PL 

x — L 


dx 2 E * Ia * L 
d 2 v p 3 * L 2 


(‘-Sr*) 

x — L 


dx 2 E * Ia (J3L — (1 — a)x ) 3 
Resolviendo por integraciones sucesivas se tiene que: 


dv p 3 * L 2 


dx E * Ia 


2 (a — l)x + (/? — a + 1 )L 
2{a — l) 2 ((a — l)x + pi) 2 


+ A 


P 3 *l 2 

v(x) = -* 


((2a — 2)x + 2 PL) ln((a — l)x + /?i) + (/? + a — 1 )L 


E * Ia 

+ Ax + B 

Para: pL<x < L 


2 (a — 1 ) 3 ((a — l)x + /?i) 


o x 

d v L ~ 1 


dx 2 E * a 3 * Ia 
d 2 v 1 

——r = ---(x - L) 

dx z E * a 6 * Ia * L 

Calculando por integraciones sucesivas se tiene que: 


dv 


dx E * a 3 * Ia * L 

1 


— Lx + C 


/ y .3 % 2 

v(x ) =---(- L - 1 - Cx + D 

E * a 3 * Ia * L\6 2 

Para encontrar las constantes de integración se toma condiciones de 


frontera. 
x = 0; v = 0 


(2 PL) ln (PL) + (p + a - 1)L 

0 —--^-1" B 

2(a — l) 3 (J3L) 

(2 PL) ln (PL) + (p + a - 1)L 


B = 


2(a — l) 3 (PL) 


x = L; v = 0 


Sílaa 3 


+ CL + D = 0 
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x = (3L;61 = 92 


91 = 


E * Ia 


2(a — l)jS + (J3 — a + 1) 
2 (a — l) 2 (a) 2 


+ A 


92 - 


E * a 6 * Ia * L \ 2 


pi 2 o 
— ~Pl 2 + C 


x = (3L;v 1 = v2 


L 3 p 3 L 3 /? 2 

= -6- + CL/3 + D 

E\aLa á 


v2 

2 3 |n (/.(a-1)/? + /./?)_ ip+Lu-L _ 

p 1 a 3 - 3a 2 + 3a - 1 ^ L(2a 4 - 8a 3 + 12a 2 - 8a + 2)/? + (2La 3 - 6La 2 + 6La - 2L)/T 

~ Êü 

+ ALP + B 


Resolviendo el sistema de ecuaciones lineales, se obtiene las siguientes 
constantes de integración: 

A = —(/? 3 * (L * a 3 * (11 — 6 * Zn(L * (a — 1) * p + L * /?)) — 18 * L * a 2 + 9 

* L * a — 2*L) + 6*L*a 3 */? 3 * Zn(L * /?) + (3 * L * a 4 — 21 

* L * a 3 + 36 * L * a 2 — 24 * L * a + 6 * L) * /? 2 + (9 * L * a 3 

— 24*L*a 2 + 21*L*a — 6 * L) * p — 2*L* a 3 + 6*L* a 2 

— 6 * L * a + 2 * L)/(6 * E * Ia * a 6 — 18 * E * Ia * a 5 + 18 * £ 

* Ia * a 4 — 6 * E * Ia * a 3 ) 


B = (2 * L 2 * p 3 * ln(L * P) + L 2 * p 3 + (L 2 * a — L 2 ) * /? 2 )/( 2 * E * Ia* a 3 
— 6 * E * Ia * a 2 + 6 * E * Ia* a — 2 * E * Ia ) 


C = —{JP 3 * (L * a 3 * (11 — 6 * Zn(L * (a — 1) * /? + L * /?)) — 18 * L * a 2 + 9 

* L * a — 2*L) + 6*L*a 3 *p 3 * ln(L * /?) + (3 * L * a 4 — 12 
*L*a 3 +18*L*a 2 — 12*L*a4-3*L)*/? 2 — 2 * L * a 3 

+ 6 * L * a 2 — 6 * L * a + 2 * L)/(6 * E * Ia * a 6 — 18 * E * Ia 

* a s + 18 * E * Ia * a 4 — 6 * E * Ia * a 3 ) 

D = (J3 3 * (i 2 * a 3 * (11 — 6 * Zn(L * (a — 1) * /? 4- L * /?)) — 18 * L 2 * a 2 + 9 

* L 2 * a — 2* L 2 ) + 6 * L 2 * a 3 * p 3 * ln(L * /?) + (3 * L 2 * a 4 
— 12 * L 2 * a 3 + 18 * L 2 * a 2 — 12 * L 2 * a + 3 * L 2 ) * p 2 )/(6 
*E * Ia * a 6 — 18 *E*Ia*a s + 18 *E*Ia* a 4 — 6*E*Ia 

* a 3 ) 


Por lo tanto, las funciones de forma resultantes son: 
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• Para: 0 < x < pL 


01 oo = 


/? 3 * L 2 
E * Ia 


((2 a — 2)x + 2/?i) ln((a — l)x + /3L ) + (/? + a — 1 )L 
2(a — l) 3 ((a — l)x + /?l) 


+ Ax + B 


• Para: (SL< x < L 


(17) 


01 oo 


I 

E * a 3 * Ia * L 



+ Cx + D 


(18) 


De esta forma las ecuaciones que gobiernan los giros de la viga son: 

• Para: 0 <x < fíL 


d01 _p 3 *L 2 
dx E * Ia 

• Para: fiL<x<L 


2 (a — l)x + (/? — a + 1 )L 
2(a — l) 2 ((a — l)x + /?l) 


d01 

dx 


1 

E * a 3 * Ia* L 



+ C 


(19) 


( 20 ) 


3.2 Funciones de forma para carteia al lado derecho 

Considerando que la viga no es simétrica, es necesario encontrar la 
segunda funciónde forma, para ello se aplica un momento unitário en ei nudo final 
de la viga, como se indica en la figura 6. 


í tf,L "'\l 

•Wb. 
l/l I 

Figura 6 Aplicación de un momento unitário en ei nudo final de una viga recta-con 
carteia a la derecha. 

M(x) = X - 
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• Para: 0 < x < pL 

, X 

d V _ _L_ 

dx2 p j (. 1 — a \ 3 

d 2 v 1 x 

dx 2 E * Ia * L / 1 — a \ 3 

d 2 y /? 3 * L 2 x 

dx 2 E * Ia (J3L — (1 — a) * x) 3 

Calculando por integrales sucesivas: 

dv 
dx 
= A 

_ L 2 */? 3 *(L*p + x* (2* a-2)) _ 

— E * Ia* ((2 * L 2 * a 2 — 4 * L 2 * a + 2 * L 2 ) * p 2 + x * (4 * L * a 3 — 12 * L * a 2 + 12 * L * a — 4 * L) * (3 + x 2 * (2 * a 4 — 8 * a 3 + 12 * a 2 — 8 * a + 2)) 


,2 , n, , / fa(L*/? + X*(«-l)) _W?_ 

^ ya 3 — 3*cr 2 + 3*a: — 1 (2* L* a 2 — 6* L* a 2 + 6*L*a — 2*L)*/?+x*(2*a 4 — 8 * a 3 + 12 * a 2 — 8 * cr + 2) 

E * Ia 

+A * x+B 

• Para: pL<x <L 

9 * 

dv = L 

dx 2 E * a 3 * Ia 


d 2 v 
dx 2 


1 

E * a 3 * Ia * L 


00 


Resolviendo por integraciones sucesivas se tiene que: 

dv x 2 

- = - -2 + C 

dx 2*E*la*L*a 


v(x) = 


x 


3 


6 * E * Ia * L * a 3 


+ Cx + D 


Para encontrar las constantes de integración se toma condiciones de 

frontera. 


x = 0 ; v = 0 
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l 2 p 3 


B — 


ln m 
3 a 2 + 3a 


l + 2LÍp- 


6 La 2 + 6 La — 2 L 


E * Ia 


= 0 


x — L;v — 0 


L 2 

6Elaa 3 


+ CL +D = 0 


x = pL)81 = 92 

(2 La — L)/? 2 — 2AE\aa 4 + AAEXaa 3 — 2AElaa 2 L/3 2 + 2CElaa 3 
2E\aa 4 — AE\aa 3 + 2£'Iacr 2 2E\aa 3 


x = pL;v 1 = v2 


L 2 /? 3 + BCElaLa 3 + 6DElaa 3 
6E\aa 3 

2 L 2 a(3 3 ln (La/?) + L 2 /? 3 + ( -2AElaLa 4 + ÒAElaLa 3 - BAE\aLa 2 + 2AE\aLd)P - 2BE\aa i + ÒBElaa 3 - ÒBElaa 2 + 2BE\aa 

2E\aa i — 6i?Iaa 3 + 6i?Iaa 2 — 2E\aa 


Calculando el sistema de ecuaciones lineales, se tiene que: 

A 

_ 6La 3 p 3 ln(LaP) + /? 3 (L« 3 (-61n(L/J) - 11)+ 18La 2 - 9La + 2L) + [9La 3 - 18La 2 + 12La — 3L)/? 2 — La 3 + 3La 2 — 3La + L 
~ 6Elaa 6 — 18£Iaa 5 + 18£Iaa 4 - 6£Iaa 3 


L 2 /? 3 (21n(L/?) + 1) 

2£laa 3 — 6£laa 2 + 6£laa — 2£la 

C 

6La 3 /? 3 ln(La/?) + /? 3 (La 3 (—6 ln(L/?) — 11) + 18La 2 — 9 La + 2 L) — La 3 + 3 La 2 — 3 La + L 
òElaa 6 — 18£laa 5 + 18£laa 4 — 6£laa 3 


D 

6L 2 a 3 /? 3 ln(La/?) +/? 3 (L 2 a 3 (-6 ln(L/?) - 11)+ 18L 2 a 2 - 9L 2 a + 2L 2 ) 
6£laa 6 — 18£laa 5 + 18£laa 4 — 6£laa 3 


Por lo tanto, las funciones de forma resultantes son: 

• Para: 0 < x < flL 


02(x) 

1 ? , n, , ( ln{L*p + x*{a- 1)) _W?_ 

^ ya 3 — 3*a 2 + 3*a — 1 (.2 * L* a 3 — 6* L* a 2 + 6*L*a — 2*L)*(3+x*(2*a 4 — 8 * a 3 + 12* a 2 — 8 * a + 2) 

E * Ia 

+ A* x + B 


Para: (SL<x<L 


( 21 ) 
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02 O) = 


+ Cx + D 


6*E* la *L*a 3 

De esta forma las ecuaciones que gobiernan los giros de la viga son: 

• Para: 0 < x < fiL 


( 22 ) 


ÍÍ01 


dx 
= A 

L 2 * fi 3 * (i * fi +x * (2 * a -2)) 

E * Ia* ((2 *L 2 *a 2 -4*L 2 *a+2*L 2 )*/? 2 +x*(4*L*a 3 -12*L*a 2 + 12*L*a-4*L)*/?+;t 2 *(2*a 4 -8*a: 3 + 12*a 2 -8*a + 2)) 


• Para: (3L< x < L 


(23) 


d01 x 2 

dx 2*E*la*L*a 3 


(24) 


3.3 Tablas para calcular matriz de flexibilidad y rigidez en sistema 1 

La tabla 3 sirve para el cálculo de los coeficientes a*,a'*,£*, con los que se 
hallan los elementos de la matriz de flexibilidad de una viga con carteia a la izquierda 
y la misma tabla sirve para calcular los coeficientes a'*,a*,£* en una viga con carteia 
a la derecha. Para cada fila hay tres valores y corresponden a los indicados. 


/ll_ EL 


fl2 — a 


EL 


(25) 


fl2 — /21 — £ p, 

t l 

Ahora, en la Tabla 4 se presentan los coeficientes: k*,k'*,a*, para el cálculo 
de los elementos de la matriz de rigidez en sistema 1, de una viga con carteia a la 
izquierda; para el caso de carteia a la derecha el orden de los coeficientes cambia 
a: k'*,k*,a* 
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Tabla 3 Coeficientes a*,a'*,E*, para ei caso de carteia ai lado izquierdo y 
coeficientes a'*,a*,E*, para ei caso de carteia ai lado derecho. 


p\a 

0,95 

0,9 

0,85 

0,8 

0,75 

0,7 

0,65 

0,6 

0,55 

0,5 

0,45 

0,4 


0,3074 

0,2826 

0,2591 

0,2368 

0,2157 

0,1958 

0,1771 

0,1597 

0,1434 

0,1283 

0,1144 

0,1017 

0,5 

0,3318 

0,3301 

0,3285 

0,3268 

0,3251 

0,3233 

0,3215 

0,3197 

0,3178 

0,3158 

0,3138 

0,3117 


-0,1620 

-0,1574 

-0,1528 

-0,1482 

-0,1437 

-0,1392 

-0,1347 

-0,1303 

-0,1260 

-0,1217 

-0,1175 

-0,1133 


0,3091 

0,2860 

0,2640 

0,2431 

0,2233 

0,2045 

0,1869 

0,1704 

0,1549 

0,1406 

0,1273 

0,1151 

0,45 

0,3322 

0,3310 

0,3298 

0,3286 

0,3273 

0,3260 

0,3247 

0,3234 

0,3220 

0,3206 

0,3191 

0,3176 


-0,1628 

-0,1589 

-0,1550 

-0,1512 

-0,1474 

-0,1436 

-0,1398 

-0,1361 

-0,1325 

-0,1288 

-0,1252 

-0,1217 


0,3110 

0,2897 

0,2694 

0,2500 

0,2316 

0,2142 

0,1978 

0,1823 

0,1679 

0,1544 

0,1418 

0,1303 

0,4 

0,3325 

0,3317 

0,3309 

0,3300 

0,3291 

0,3282 

0,3273 

0,3263 

0,3254 

0,3244 

0,3233 

0,3223 


-0,1635 

-0,1603 

-0,1571 

-0,1540 

-0,1509 

-0,1478 

-0,1447 

-0,1417 

-0,1386 

-0,1356 

-0,1327 

-0,1297 


0,3131 

0,2937 

0,2752 

0,2576 

0,2408 

0,2249 

0,2098 

0,1956 

0,1823 


0,1582 


0,35 

0,3328 

0,3322 

0,3317 

0,3311 

0,3305 

0,3299 

0,3293 

0,3286 

0,3280 


0,3266 



-0,1641 

-0,1616 

-0,1591 

-0,1567 

-0,1542 

-0,1517 

-0,1493 

-0,1469 

-0,1444 

-0,1421 

-0,1397 



0,3154 

0,2981 

0,2817 

0,2659 

0,2509 

0,2367 

0,2232 

0,2104 

0,1983 


0,1765 

0,1667 

0,3 

0,3330 

0,3326 

0,3323 

0,3319 

0,3316 

0,3312 

0,3308 

0,3304 

0,3300 


0,3291 

0,3287 


-0,1648 

-0,1629 

-0,1610 

-0,1591 

-0,1572 

-0,1553 

-0,1535 

-0,1516 

-0,1498 

-0,1480 

-0,1461 

-0,1443 


0,3178 

0,3029 

0,2887 

0,2750 

0,2620 

0,2496 

0,2378 

0,2266 

0,2161 


0,1968 

0,1881 

0,25 

0,3331 

0,3329 

0,3327 

0,3325 

0,3323 

0,3321 

0,3319 

0,3316 

0,3314 


0,3309 

0,3306 


-0,1653 

-0,1639 

-0,1626 

-0,1612 

-0,1599 

-0,1585 

-0,1572 

-0,1559 

-0,1545 

-0,1532 

-0,1519 

-0,1506 


0,3205 

0,3081 

0,2963 

0,2849 

0,2741 

0,2637 

0,2538 

0,2445 

0,2356 

0,2272 

0,2193 

0,2119 

0,2 

0,3332 

0,3331 

0,3330 

0,3329 

0,3328 

0,3327 

0,3326 

0,3325 

0,3323 

0,3322 

0,3321 

0,3319 


-0,1658 

-0,1649 

-0,1640 

-0,1631 

-0,1622 

-0,1613 

-0,1604 

-0,1595 

-0,1587 

-0,1578 

-0,1569 

-0,1561 


0,3234 

0,3137 

0,3045 

0,2957 

0,2872 

0,2791 

0,2713 

0,2640 

0,2570 

0,2504 

0,2441 

0,2382 

0,15 

0,3333 

0,3332 

0,3332 

0,3332 

0,3331 

0,3331 

0,3330 

0,3330 

0,3329 

0,3329 

0,3328 

0,3327 


-0,1661 

-0,1656 

-0,1651 

-0,1646 

-0,1641 

-0,1636 

-0,1630 


-0,1620 

-0,1615 

-0,1610 

-0,1605 


Los elementos de la matriz de rigidez, vienen dados por las siguientes 
ecuaciones: 


fcn = k * 


EI, 


o 


k 22 = k' 


L 

El n 


(16) 


~ ^21 — a 


EI n 


En la tabla 4, se presentan los coeficientes que permiten calcular los 
elementos de la matriz de rigidez cuando se tiene carteia al lado derecho de la viga. 
































24 


Roberto Aguiar y Alejandro Bonilla 


Tabla 4 Coeficientes k*,k'*,a*, para el caso de carteia al lado izquierdo y 
coeficientes k'*,k*,a *, para el caso de carteia al lado derecho. 


p\a 

0,95 

0,9 

0,85 

0,8 

0,75 

0,7 

0,65 

0,6 

0,55 

0,5 

0,45 

0,4 


4,3806 

4,8160 

5,3165 

5,8954 

6,5687 

7,3567 

8,2852 

9,3864 

10,7014 

12,2820 

14,1936 

16,5171 

0,5 

4,0588 

4,1231 

4,1937 

4,2719 

4,3587 

4,4557 

4,5649 

4,6884 

4,8290 

4,9901 

5,1754 

5,3896 


2,1390 

2,2954 

2,4722 

2,6731 

2,9026 

3,1665 

3,4717 

3,8268 

4,2427 

4,7326 

5,3128 

6,0028 


4,3596 

4,7674 

5,2316 

5,7625 

6,3726 

7,0766 

7,8931 

8,8441 

9,9563 

11,2620 

12,7985 

14,6085 

0,45 

4,0569 

4,1192 

4,1876 

4,2631 

4,3466 

4,4395 

4,5433 

4,6596 

4,7905 

4,9382 

5,1050 

5,2934 


2,1360 

2,2881 

2,4588 

2,6511 

2,8689 

3,1164 

3,3990 

3,7230 

4,0956 

4,5256 

5,0226 

5,5972 


4,3354 

4,7118 

5,1355 

5,6140 

6,1561 

6,7720 

7,4738 

8,2751 

9,1919 

10,2416 

11,4429 

12,8146 

0,4 

4,0551 

4,1152 

4,1809 

4,2530 

4,3324 

4,4199 

4,5166 

4,6237 

4,7423 

4,8739 

5,0194 

5,1799 


2,1314 

2,2771 

2,4391 

2,6200 

2,8223 

3,0493 

3,3046 

3,5922 

3,9166 

4,2827 

4,6953 

5,1588 


4,3077 

4,6489 

5,0281 

5,4502 

5,9208 

6,4464 

7,0338 

7,6907 

8,4249 

9,2443 

10,1564 

11,1669 

0,35 

4,0529 

4,1102 

4,1725 

4,2403 

4,3141 

4,3947 

4,4825 

4,5783 

4,6827 

4,7962 

4,9189 

5,0510 


2,1247 

2,2617 

2,4125 

2,5786 

2,7620 

2,9647 

3,1888 

3,4365 

3,7102 

4,0119 

4,3433 

4,7054 


4,2763 

4,5785 

4,9093 

5,2718 

5,6690 

6,1042 

6,5808 

7,1021 

7,6711 

8,2902 

8,9607 

9,6827 

0,3 

4,0499 

4,1036 

4,1614 

4,2235 

4,2904 

4,3624 

4,4397 

4,5225 

4,6111 

4,7054 

4,8052 

4,9101 


2,1158 

2,2415 

2,3781 

2,5265 

2,6878 

2,8630 

3,0532 

3,2593 

3,4820 

3,7219 

3,9789 

4,2523 


4,2409 

4,5003 

4,7796 

5,0802 

5,4035 

5,7507 

6,1228 

6,5206 

6,9446 

7,3944 

7,8693 

8,3670 

0,25 

4,0459 

4,0948 

4,1467 

4,2018 

4,2602 

4,3220 

4,3873 

4,4560 

4,5281 

4,6032 

4,6810 

4,7609 


2,1043 

2,2160 

2,3355 

2,4633 

2,5998 

2,7455 

2,9004 

3,0649 

3,2388 

3,4217 

3,6130 

3,8116 


4,2013 

4,4144 

4,6396 

4,8774 

5,1279 

5,3913 

5,6675 

5,9560 

6,2563 

6,5674 

6,8877 

7,2155 

0,2 

4,0406 

4,0831 

4,1276 

4,1742 

4,2227 

4,2732 

4,3255 

4,3795 

4,4350 

4,4917 

4,5492 

4,6071 


2,0900 

2,1847 

2,2845 

2,3893 

2,4991 

2,6140 

2,7338 

2,8582 

2,9869 

3,1192 

3,2546 

3,3920 


4,1575 

4,3210 

4,4904 

4,6654 

4,8460 

5,0317 

5,2220 

5,4164 

5,6142 

5,8144 

6,0159 

6,2176 

0,15 

4,0336 

4,0682 

4,1038 

4,1403 

4,1777 

4,2159 

4,2547 

4,2940 

4,3336 

4,3733 

4,4128 

4,4518 


2,0725 

2,1476 

2,2251 

2,3049 

2,3869 

2,4710 

2,5568 

2,6440 

2,7323 

2,8213 

2,9104 

2,9990 


4. NUEVAS FUNCIONES DE CEINCI-LAB 

Se presentan los programas con los cuales se halla en primer lugar la matriz 
de rigidez en coordenadas globales, para una viga con carteias, y el programa que 
encuentra la contribución de las vigas con carteias a la matriz de rigidez de la 
estructura, todo esto para una estructura en que no se considera nudo en la mitad 
de la viga. 

• [K3]=kmiembro_cartelalineal(b,h,H,L,Lv,E,seno,coseno) 

Los datos son: b es el ancho constante de la viga; h es la altura de la viga 
en el centro de luz; Hes la altura de la viga en el nudo inicial; L es la luz de 
la viga; Lv es la longitud de la carteia; £ es el módulo de elasticidad dei 
material; seno es el seno dei ângulo que forma el eje de la viga con el eje 
X; coseno es el valor dei coseno. 

• [SS]=krigidez_cartela(ngl,EL EM, L, se no, coseno, VC,E) 
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Los datos son: ngl número de grados de libertad; ELEM es una matriz que 
contiene para cada viga con carteias: b, h, H, Lvque son la base, altura en 
centro de luz; altura en nudo inicial y longitud de la carteia. Son 4 datos que 
tiene ELEM para cada elemento; L, seno, coseno, son vectores con la 
longitud, seno y coseno de cada elemento; VC matriz que contiene los 
vectores de colocación de cada elemento; E es el módulo de elasticidad dei 
material. 

Ahora para el caso en que se considere nudo en la mitad de la viga, los 
programas que encuentran la matriz de rigidez de un elemento y la contribución a 
la matriz de rigidez de la estructura, son: 

• [K3]=kmiembro_cartelalineal_nudointermedio(b,h,H,L,Lv,E,seno,cose 
no,tipo) 

Los datos son: b es el ancho constante de la viga; h es la altura de la viga 
en el centro de luz; H es la altura de la viga en el nudo inicial; L es la luz 
para el elemento cartela-viga de sección constante, es decir vale la luz de 
la viga para 2; Lves la longitud de la carteia; E es el módulo de elasticidad 
dei material; seno es el seno dei ângulo que forma el eje de la viga con el 
eje X; coseno es el valor dei coseno; tipo factor que identifica si la carteia 
se halla a la izquierda o a la derecha, para el primer caso tipo=0, para el 
segundo caso tipo=1. 

• [SS]=krigidez_cartela_nudointermedio(ngl, ELEM, L, seno, cose, VC,E) 

Los datos son: ngl número de grados de libertad; ELEM es una matriz que 
contiene para cada viga con carteias: b, h, H, Lv,tipo que son la base, altura 
en centro de luz; altura en nudo inicial, longitud de la carteia y el tipo que es 
igual a 0 si la carteia se halla a la izquierda e igual a 1 si la carteia se 
encuentra a la derecha. Son 5 datos que tiene ELEM para cada elemento; 
L, seno, cose, son vectores con la longitud dei elemento hasta la mitad de 
la viga, seno y coseno de cada elemento; VC matriz que contiene los 
vectores de colocación de cada elemento; E es el módulo de elasticidad dei 
material. 

. [ELEMVC]=gelem_portico_cartela(SECCION,UBICACION) 

El programa gelem_portico_cartela genera el archivo ELEMVC que 
contiene la siguiente información, para cada elemento de una viga, con 
carteia (de preferencia) o sin carteia: 

b h H Lv Tipo 

Se genera para cada elemento, se recuerda que Tipo=0 para el caso de que 
carteia se halle a la izquierda y vale 1 cuando carteia se encuentra a la 
derecha. Los datos de entrada son dos matrices denominadas SECCION, 
UBICACION. 
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Para la matriz SECCION los datos son: b ancho de la viga con carteia; h 
altura en el centro de luz de la viga; H altura en el extremo de la viga con 
carteia; Lv es la longitud de la carteia; la primera fila de datos es para la 
sección 1, la segimda fila para la dos, etc., esta información se da en la 
matriz UBICACIÓN que tiene 5 datos por cada fila y son los siguientes: 

i) Identificación de la sección de la viga. 

ii) Identificación dei elemento 

iii) Número de elementos a generar (que tienen la misma sección) 

iv) Incremento en la numeración de los elementos 

v) El código para saber si carteia está a la izquierda (0) o a la derecha 

( 1 ) 


Los programas dei sistema de computación CEINCI-LABse están subiendo 
a la partición dei primero de los autores en Research Gate, para estructuras con 
carteias se hallan en la carpeta: Análisis Plano Básico. También se hallan 
publicadas en Aguiar et ai (2019, 1,2). 


5. EJEMPLO DE APLICACIÓN 

La estructura de la figura 6 a,b y c; es simétrica en planta y elevación, con 
luces de 5 y 10 metros, y con vigas rectangulares de 40/60 para las luces de 5.0 m, 
y vigas con carteias lineales para las luces de 10.0 m. La altura de cada piso es de 
3.60 m. La carga muerta D = 0.8 T/m 2 , en los pisos 1 al 3 y D = 0.6 T/m 2 , en el piso 
4. La carga viva L = 0.2 T/m 2 , en los tres primeros pisos y la cubierta es inaccesible. 
Se pide: 


i) Encontrar la matriz de rigidez lateral de un pórtico tipo exterior, calcular 
con un módulo de elasticidad E h = 2400000 T/m 2 (Figura 7) 

ii) Encontrar la matriz de rigidez en coordenadas de pisos, considerando 
un grado de libertad por piso, ubicado en el Centro de Masas, de cada 
planta. (Figura 6 b) 

iii) Flallar la matriz de Masas para el modelo de un grado de libertad por 
piso, teniendo en cuenta la presencia de la grada en los 3 primeros 
pisos. (Figura 6 c) 

iv) Resolver el problema de valores y vectores propios. 

v) Encontrar un modelo equivalente de un grado de libertad. 

vi) Realizar un análisis lineal y no lineal en el sistema equivalente de un 
grado de libertad, ante la componente N-S, dei registro de 
aceleraciones obtenido en Manta, dei terremoto dei 2016. Considerar, 
que el desplazamiento de fluência q y = 0.001 * h eq . Donde h eq es la 
altura equivalente dei sistema de un grado de libertad. 

vii) Encontrar la respuesta de desplazamientos en cada uno de los pisos 
de la estructura. 
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(b) 



(c) 

Figura 6 Estructura de análisis; a) Vista en 3D; b) Modelo espacial de 1 grado de 
libertad por planta; c) Vista en planta. 
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Figura 7 Geometria de Pórtico Exterior Tipo, columnas exteriores de 60/60 cm y 
columnas interiores de 110/110 cm. 

• SOLUCIÓN 


i) Cálculo de matriz de rigidez lateral de Pórtico Exterior 


Lo más fácil es suponer que las vigas rectangulares de 40/60 cm, indicadas 
en la figura 7, son vigas con carteia, para esto se tiene: 

h = H = 0.60m b — 0A0m ; L v =0.0m 

De esta manera la matriz de rigidez de todo ei pórtico que se denomina KT, 
es igual a la matriz de rigidez de las columnas rectangulares KH, más la matriz de 
rigidez de las vigas rectangulares y con carteias (calculadas todas como con carteia) 

KVC. 


Por otra parte, se va a trabajar con un modelo que considera nudo en la 
mitad de la viga. El usuário da los datos de la separación de los vanos sv , y de la 
separación de las alturas de pisos sp, con estos dos vectores ya se puede dibujar 
la estructura que apare en la figura 8 a, con la numeración de nudos y elementos. 
Esta gráfica es importante imprimiria para la entrada de los datos de las secciones 
de las columnas y vigas; se aprecia que dei elemento 1 ai 16 son columnas y dei 17 
ai 40 son vigas. Para facilitar la entrada de datos de las vigas con carteia, se tiene 
el programa denominado gelem_portico_cartela. 
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En la figura 8 b, se muestran los grados de libertad dei modelo de cálculo, 
con azul se presentan las coordenadas principales o coordenadas laterales para las 
cuales se condensa la matriz de rigidez. 


Numeración de nudos y elementos 



(a) 


Grados de Libertad 



(b) 

Figura 8 a) Numeración de nudos y elementos, hallado con la función dibujoplano, 
b) Numeración de los grados de libertad, con azul se indican las 
coordenadas principales, encontrado con la función dibujogdl 
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El listado dei programa que conduce al cálculo de la matriz de rigidez lateral 
de un pórtico interior se indica a continuación (utilizando las funciones de CEINCI- 
LAB) 


% Ej_7_clase_ESPE 
clear all; cic 
%% Pórtico 1 y 4 

sv=[5.0;10.0;5.0]; % luz dei vano en metros 
sp=[3.6;3.6;3.6;3.6]; % altura de cada piso en metros. 
nod=32; nr=4; E=2400000; 

Seccion=[1 0.60.6 3 4; 2 1.10 1.10 34; 

3 1.1 1.1 3 4;4 0.6 0.6 3 4]; 

%% PÓRTICO SIN DIAGONALES 


[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp); 
[X,Y]=glineajoortico2(nv,np,sv,sp,nod,nr); %coordenadas 
[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 

%dibujoplano(X,Y,NI,NJ) 

%% Geometria de la estructura 
[L,seno,coseno]=longitud (X,Y,NI,NJ); 

[CG,ngl]=cg_sismo2(nod,nr,Y); 

%dibujogdl(X,Y,NI,NJ,CG); 

[VC]=vc(NI,NJ,CG); 

% Contribución de columnas rectangulares 

LC=L( 1:16) ;senC=seno( 1:16) ;cosC=coseno( 1:16); VCC=VC( 1:16,:); 

[ELEMC]=gelem_portico(Seccion); % tiene base y altura 

[KC]=krigidez(ngl,ELEMC,LC,senC,cosC,VCC,E); 

% Contribución de vigas rectangulares y con carteias 
b1=0.40;h1=0.60;H1=0.60;LV1=0.0;% Son datos de viga rectangular 
b2=0.40;h2=0.60;H2=1,00;LV2=2.00;% Viga con carteia 
SECCIONVC=[b1 hl Hl LV1;%[b h H LV] 
b2 h2 H2 LV2]; 

UBICACION=[1 17 3 6 0; %[tipo de seccion, elemento, #generar,inc, lado carteia] 


1 1836 1 


% lado=0 izquierda; lado=1 derecha 


2 19360; 

2 20 3 6 1; 

1 21 3 6 0; 

1 22 3 6 1]; 

[ELEMVC]=gelem_portico_cartela(SECCIONVC, UBICACION); 

LVC=L( 17:40) ;senVC=seno( 17:40) ;cosVC=coseno( 17:40); VCVC=VC( 17:40,:); 
[KVC]=krigidez_cartela_nudointermedio(ngl,ELEMVC,LVC,senVC,cosVC,VCVC,E) 


KT=KC+KVC; 

%% Condensación estática de K 
na=4; 

kaa=KT (1: na, 1: na); kab=KT( 1: na, na+1: ngl); kba=kab'; 
kbb=KT( na+1: ng I, na+1: ngl); 

KL=kaa- kab* i nv( kbb) * kba 


■ 211590 -129170 

„(i) _ «.(4) _ -129170 167500 

~ n L ~ 44320 -112180 

-6365.6 29642 


44320 -6365.6- 

-112180 29642 

125980 -49329 

-49329 24784 - 
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La matriz de rigidez lateral de un Pórtico interior es diferente ya que todas 
las columnas son de 1.10/1.10 m. 


■ 382400 
k ( 2 ) _ „( 3 ) _ -233880 
i i 80910 

- -11540 


-233880 

301720 

-202500 

53860 


80910 -11540- 

-202500 53860 

223610 -86030 

-86030 41420 - 


ii) Matriz de rigidez en coordenadas de piso 


Para ei modelo de un grado de libertad por piso, la matriz de rigidez, que se 
denomina K xx es igual a la suma de las matrices de rigidez lateral de los pórticos 
en sentido X. 


K xx = k[ 1] 

+ K ( L l) + K 

?+*í 4) 


- 1.1880 

-0.7261 

0.2505 

-0.0358 

-0.7261 

0.9384 

-0.6294 

0.1670 

0.2505 

-0.6294 

0.6992 

-0.2707 

— 0.0358 

0.1670 

-0.2707 

0.1324 


iii) Matriz de masas en coordenadas de piso 


-33.31 0 0 0 

0 33.31 0 0 

0 0 33.31 0 

0 0 0 24.49- 


%% Matriz de Masas, Períodos y Modos de vibración 
D=0.8;L=0.2;Pt=D+0.25*L 
area=20*20-4*4; %Calculo sin grada 
W1=area*Pt;W2=W1 ;W3=W1; 

Pt4=0.6; %Losa inaccesible 
W4=20*20*Pt4; 

m1=W1/9.8;m2=W2/9.8;m3=W3/9.8;m4=W4/9.8; 

M=mdiag(m1,m2,m3,m4) 

[T,fi,OM]=orden_eig(Kxx,M) 


iv) Períodos y Modos de vibración 


En la tabla 5 se indican los períodos de vibración y en la figura 9 los modos 
de vibración encontrados con ei programa orden_eig de CEINCI-LAB 


Tabla 5 Períodos de vibración encontrados con programa orden_eig 


Modo 

1 

2 

3 

4 

T (s) 

0.4790 

0.1031 

0.0415 

0.0249 















32 


Roberto Aguiar y Alejandro Bonilla 


0.017' 


' 0.066 ' 


' 0.117 ' 


' 0.108 ' 

0.056 

0 ( 2 ) = 

0.118 

0 ® = 

0.027 

0« = 

- 0.111 

0.102 

0.046 

-0.109 

0.074 

0.148. 


-0.115. 


. 0.070 . 


-0.030. 



0030 

M - 



Cuarto modo 


Figura 9 Modos de vibración de modelo de un grado de libertad por piso en X 


v) Modelo equivalente de un grado de libertad 

Para ei cálculo dei modelo de un grado de libertad, se trabaja solo con el 
primer modo. Luego T 1 = 0.479 s. Las ecuaciones de cálculo dei sistema 
equivalente de un grado de libertad, al sistema de 4 grados de libertad, se encuentra 
con las siguientes ecuaciones. Aguiar (2012) 

_ (0 (1)£ M [l]) 2 

meq ~ M 0W 



c eq = 2Ç \m eQ k , 


eq ,x -eq 


(0 (1)t M/i) 


heq 0(Dt M [1] 
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Donde M es la matriz de masas dei sistema de n grados de libertad; T es el 
período fundamental dei sistema de n grados de libertad que se halla resolviendo el 
problema de valores y vectores propios con My K. Siendo K la matriz de rigidez dei 
sistema de n grados de libertad; [1] es un vector unitário que tiene n unos, para la 
estructura de la figura 6, tendrá 4 unos. (f> m es el primer modo de vibración de la 
estructura con n grados de libertad; h es un vector con las alturas a cada piso, 
medido desde la base. 

Al utilizar las ecuaciones de cálculo para hallar el sistema equivalente de un 
grado de libertad, se obtiene: 


m eq = 89.617 T.s 2 /m 
k eq = 15418 T/m 
c eq = 11.039 T.s/m 
H eq = 9.47 m 
f\ = 9.4667 

Las sentencias con las cuales se hallan el sistema equivalente de un grado 
de libertad se indican a continuación. 

%% Modelo equivalente de 1 gdl 
T1=T(1);fi1=fi(:,1); 

uno=[1 ;1 ;1 ;1];hh=[3.6;2*3.6;3*3.6;4*3.6] 

meq=(f i 1 ’*M*uno) A 2/(fi 1 '*M*fi 1) 

keq=meq*(2*pi/T1) A 2 

ceq=2*0.05*sqrt(meq*keq) 

heq=(fi1'*M*hh)/(fi1'*M*uno) 

gamai =(f i 1 '*M*uno)/(fi 1 '*M*fi 1) 


vi) Análisis Lineal y no Lineal dei sistema equivalente de un grado de 
libertad 


Tanto para el análisis lineal, como no lineal se trabajó con el Método p de 
Newmark, claro está, en el análisis lineal la rigidez dei sistema se mantiene 
constante, en cambio en el análisis no lineal, esta va cambiando en función de la 
rama en que se encuentra en la constitutiva indicada en la figura 10, que fue con la 
que se trabajó (Elasto-Plasto) 

%% Datos de entrada 

load rec_manta % Registro de Manta 

% comp: v =componente vertical; ns =componente norte-sur; ew =componente 

este oeste 

comp=ns; 

M=89.617; %Ts2/m 
K=15418; %T/m 
heq=9.47; %m 
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zeda=0.05; beta=1/6; 

uy=0.00rheq; % Desplazamiento de Fluência 
gamai =9.4667; % Factor de participacion modal 
p=0.01*comp; % registro en m/s2 

g=p'; 

Dibujo=1; 

%% Respuesta elástica 

[ttn,upp,up,u,Sdmax]=beta_newmark_W(zeda,K,M,Dt,comp); 

%% Respuestas no lineal 

[ttiempo,df,vf,acf,fuerza,dmaxf,vmaxf,acmaxf]=elastoplasto_nolineal(M,K,zeda,Dt,b 

eta,uy,g,Dibujo); 




Figura 10 Constitutiva Elasto perfectamente plástico. El desplazamiento de fluência 
q y , se considero q y = 0.001 * h eq , en ei Ejemplo. 

En la figura 11, se presenta, a la izquierda la curva constitutiva dei sistema 
de un grado de libertad, ante la componente N-S dei registro de aceleraciones 
obtenido en Manta, dei terremoto de 2016; a la derecha se indica en este orden la 
respuesta de desplazamientos, de velocidad y de aceleración; con línea azul se 
tiene la respuesta lineal y con rojo la no lineal. 

El desplazamiento lateral máximo, hallado en la respuesta lineal es 
5.80 cm; y el desplazamiento máximo, encontrado en la respuesta no lineal es 
5.26 cm. 
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Figura 11 Curva constitutiva a la izquierda y respuesta lineal (azul) y no lineal (rojo) 
en desplazamiento, velocidad y aceleración, dei sistema equivalente de 1 
gdl. 


vii) Respuesta de desplazamientos en la estructura de 4 pisos 


Se tiene la respuesta no lineal en el sistema de un grado de libertad, ahora 
para pasar al sistema de múltiples grados de libertad, se lo hace con el factor de 
participación dei primer modo r\ y para encontrar los desplazamientos en cada piso 
se utiliza el primer modo de vibración 

q = ri* s dl *0« 


q = 9.4667 * q lne 


0.017' 

0.056 

0.102 

0.148. 


Donde q ine es la respuesta en el tiempo de desplazamientos que se obtiene 
en el análisis no lineal, con el programa elastoplasto_nolineal Por cierto, a esta 
variable se ha denominado df, de esta forma se halla la respuesta en cada uno de 
los pisos, la misma que se presenta en las figura 12. 
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(c) (d) 

Figura 12 Respuesta en el tiempo, en cada piso de la estructura; a) Historia de 
desplazamientos en primer piso; b) En segundo piso; c) En tercer piso; 
d) En cuarto piso. 

Las respuestas máximas, que se hallan en cada piso, en valor absoluto, se 
indican en tabla 6 


Tabla 6 Desplazamientos Máximos en cada piso 


Piso 

i 

2 

3 

4 

Rmax (CITl) 

0.9 

2.86 

5.13 

7.28 


6. COMENTÁRIOS Y CONCLUSIONES 

Se ha presentado con bastante detalle, el cálculo de las funciones de forma 
para una viga con carteias de hormigón armado; después se ha determinado los 
elementos de la matriz de flexibilidad para el sistema de coordenadas uno, que es 
un elemento apoyado en el nudo inicial y con rodillo en el nudo final. Luego se 
obtuvo la matriz de rigidez, invirtiendo la matriz de flexibilidad. 

El cálculo se lo realizo para tres modelos, el primero de ellos no contempla 
nudo en la mitad de la viga y los otros dos modelos si lo contemplan de tal manera 
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que en uno de ellos la carteia se halla a la izquierda y en el otro la carteia se 
encuentra a la derecha. 

En base a los resultados obtenidos se han desarrollado nuevas funciones 
para determinar la contribución a la matriz de rigidez de las vigas de una estructura 
con carteias; estas funciones han sido incorporadas al sistema de computación 
CEINCI-LAB y el lector puede descargar estas y otras funciones, de la carpeta 
Pórtico Plano Básico, que el primero de los autores tiene en su partición en 
Research Gate. 

Finalmente se ha realizado el análisis sísmico espacial de una estructura 
con carteias de hormigón armado, con un modelo de un grado de libertad por planta, 
encontrando primeramente un sistema equivalente de un grado de libertad, 
encontrando la respuesta lineal y no lineal ante el registro N-S, dei acelerograma 
obtenido en Manta, en el terremoto de 2016 (M=7.8). 

Desde el punto de vista estructural, las vigas con carteia es una solución 
óptima cuando se tienen luces grandes y desde el punto de vista arquitectónico este 
tipo de vigas permite el paso de las instalaciones por la parte central de la viga. 
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